NACA TN 3874 



NATIONAL ADVISORY COMMITTEE 
FOR AERONAUTICS 


TECHNICAL NOTE 3874 


EXPERIMENTAL AND PREDICTED LATERAL- DIRECTIONAL DYNAMIC- 
RESPONSE CHARACTERISTICS OF A LARGE FLEXIBLE 35° 
SWEPT-WING AIRPLANE AT AN ALTITUDE 
OF 35,000 FEET 

By Stuart C. Brown and Euclid C. Holleman 

Ames Aeronautical Laboratory 
Moffett Field, Calif. 



Washington 

December 1956 


NATIONAL ADVISORY COMMITTEE FOR AERONAUTICS 


TECHNICAL NOTE 387 ^ 


EXPERIMENTAL AND PREDICTED LATERAL-DIRECTIONAL DYNAMIC- 
RESPONSE CHARACTERISTICS OF A LARGE FLEXIBLE 35° 
SWEPT-WING AIRPLANE AT AN ALTITUDE 
OF 35,000 FEET 

By Stuart C. Brown and Euclid C. Holleman 
SUMMARY 


The lateral-directional dynamic-response characteristics of a large 
flexible airplane are presented for flight conditions of 0.6 to 0.82 Mach 
number and an altitude of 35,000 feet. Transient responses at various 
stations on the airplane to aileron and rudder pulse-type disturbances 
have been recorded and converted into frequency-response form. A fre- 
quency range of 0.5 to 25 radians per second is presented, which includes 
the Dutch roll and rolling modes and also three structural modes. The 
measured frequency responses are compared with those obtained from trans- 
fer functions based on the lateral-directional equations of motion, which 
include first-order effects of structural deflections. Reasonably good 
agreement is obtained in the low frequency range corresponding to the 
lateral-directional modes. However, in the higher frequency range, the 
effects of the structural modes are sufficiently large so that these pre- 
dictions become inadequate even for the quantities measured near the 
center of gravity. Comparisons are also made with frequency responses 
obtained from transfer- function coefficients calculated for the rigid 
airplane to indicate the effects of flexibility. 

Calculated simplified transfer- function coefficients based on the 
lateral-directional equations of motion are compared with coefficients 
obtained by matching transient time histories with an analog computer. 
Aerodynamic derivatives extracted from measured transfer- function coef- 
ficients are also compared with estimated derivatives. In general, good 
agreement between measured and estimated values was obtained. 


INTRODUCTION 


Because of the increased importance of structural flexibility on 
aircraft stability, control, and loads, the NACA is currently flight 
testing a large flexible 35° swept-wing airplane over its entire operating 
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range of flight conditions. One phase of this program is to investigate 
characteristics of the airplane dynamic response to control-surface 
motions. Some results for the longitudinal case and a limited amount 
of data for the lateral-directional case were presented in references 1 
and 2. In a separate investigation, reference 3, some measured frequency- 
response data were also reported. The purpose of the present report is to 
present a more detailed description of the lateral-directional character- 
istics of the airplane at an altitude of 35*000 feet. 

The dynamic-response characteristics can usually he separated into 
the relatively low frequency lateral-directional modes and the higher 
frequency structural modes. The lateral-directional modes are of primary 
importance in defining stability and control characteristics of the air- 
plane as a whole, while the structural modes are of importance in the 
stability of the structure by itself or in combination with an automatic 
control system (flutter). In the present report, the experimental response 
over the low frequency range is compared with the predicted responses of 
the airplane with pseudostatic effects of flexibility included. Predicted 
responses of the rigid airplane are shown to indicate the effects of flexi- 
bility. The measured dynamic response is also examined over the higher 
frequency range to investigate effects of the first three structural modes. 

With regard to measuring and analytical techniques, the frequency 
response is a convenient way of expressing the dynamic -response character- 
istics in a standard manner which is independent of the particular input 
from the control surface. Transient responses to pilot-applied pulse 
inputs through the rudder or aileron control systems were measured and 
transformed to frequency- response form by means of the Fourier integral 
(e.g., ref. 4 ). While the frequency response represents the dynamic 
characteristics of the airplane, a sometimes more useful form for the 
analysis and synthesis of automatic control systems is the transfer func- 
tion. For the present investigation, transfer-function coefficients 
obtained from the lateral-directional three-degree-of-freedom case were 
evaluated by matching time histories on an analog computer in a manner 
similar to that described in reference 4 . Corresponding predicted 
transfer-function coefficients were calculated, with the coefficients 
modified to include flexibility effects, and are compared with the 
experimental values. 

The measured transfer-function coefficients usually must be extrapo- 
lated to different flight conditions. Since this extrapolation can be 
done more easily in terms of the aerodynamic derivatives, it is desirable 
to reduce the transfer-function coefficients to derivative form wherever 
possible. Accordingly, aerodynamic derivatives obtained from simplifica- 
tions of the principal transfer-function coefficients are compared 
with estimated derivatives . 

The data used in this report were obtained by the High-Speed Flight 
Station of the NACA as part of the cooperative program with the Ames 
Aeronautical Laboratory. 
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NOTATION 



aspect ratio 

lift coefficient 

rolling-moment coefficient 

yawing-moment coefficient 

side-force coefficient 

differential operator , — 

dt 

moment of inertia about the X axis, slug-ft 2 
effective rolling moment of inertia for flexible airplane, 

T x ■ % shc i$> slu s- ft2 

moment of inertia about the Z axis , slug-ft 2 

effective yawing moment of inertia for flexible airplane, 
I Z " < J 0 sbc n..> slug-ft 2 


T xz 

product of inertia, slug-ft 2 


dimensionless radius of 
spans 

gyration about the 

K P/ 8 a 

p 

gain of simplified -f- 

5 a 

transfer function 

K r/ &r 

I* 

gain of simplified — 

5 r 

transfer function 

K p/5 r 

gain of simplified — 

8 r 

transfer function 

M 

Mach number 


S 

wing area, sc[ ft 


V 

velocity, ft/sec 



W gross weight of airplane, lb 

b wing span, ft 


axis, wing 
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local chord of the wing measured parallel to the plane of 
symmetry 

section lift coefficient 
acceleration due to gravity, ft/sec 2 
mass of airplane, slugs 

effective mass for lateral acceleration of flexible airplane, 



normal acceleration, positive downward, gravity units 

rolling velocity, radians /sec 

dynamic pressure, lb/sq ft 

yawing velocity, radians /sec 

time, sec 

structural deflection, positive downward, in. 
vertical coordinate, positive downward, in. 
angle of attack, radians 

angle of sideslip, radians (except as otherwise noted) 

total aileron deflection, measured in a plane perpendicular 
to the hinge line, positive right aileron up, radians 
(except as otherwise noted) 

rudder deflection, measured in a plane perpendicular to the 
hinge line, positive trailing edge left, radians (except 
as otherwise noted) 

damping ratio 

damping ratio of numerator term of — transfer function 

5 a 

spanwise coordinate, in wing semispans 

parameter used in evaluating aerodynamic derivatives from 
transfer-function coefficients, defined in equation (10) 

rolling-mode time constant, sec 
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/outputN 

\input/ 


angle of bank, radians 

phase angle of output quantity minus phase angle of input 
quantity 


0) 

Wo 


OJ. 


n 


angle of yaw, radians 
frequency, radians/sec 

undamped natural frequency of numerator term of transfer 

function a 

undamped natural frequency, radians/sec 

~\ 


C 7 > 

V 

Cn |3^ 

C ^8 ^ 

C *V 


%’ 


Cy. , 
> 

°V 

Cy , 

V 

C n ., 

t 

^cp 


c n „ 


* 

+ 


C 7 } 

L r 

S' 

c z , 




derivative of coefficient with respect to subscript 


x 


2V 


S 




2V ^ Z P 


C 7 , per sec 


<lo Sb 


2VTy It 

A-ni 


C ? , per sec 




cl Sb 

o ^ 2 

C, . per sec 

% V 


n t 


Nv 


q-Sb 2 

Cr, , per sec 

2VI~ V 


2VI 


% 

2 


Zp 


c n r > P er sec 



6 


NACA TN 3874 


Nr 


— — C n , per sec 2 

^ V 


q Q S 

— - Cy , per sec 

V V 


<p 


per sec 


<lo Sb „ 2 

Cj , per sec 


Nc 


_ 2 
C,, , per sec 

I 7 n 5 


<lo S 


C v * per sec 

m„V y 6’ 


t q_SbC 7 .. 

T XZ L J/ 


L X f 


- L xz 

I z F 


V 


% + r X F B P 


N 1 

r 


L 3 + r X F N 3 


W r + r Z F L r 


V 

V 


W P + % L P 

^ + r X F N 6 


N ' 

5 


W 6 + r Z F L 5 
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Subscripts 


B body axis 

eg center of gravity 

t tail 

wt left wing tip 

Dots are used to indicate differentiation with respect to time; for 
example, p = 


TEST EQUIPMENT 


The test airplane was a Boeing B-^7A with wing vortex generators 
installed (fig. 1). An external nose boom and an optigraph for measuring 
the movements of target lights mounted on the wing and tail were added 
(figs. 1 and 2). Pertinent physical characteristics are listed in table I 
Left and right aileron, and rudder deflections were each measured by three 
NACA resistance-type control-deflection indicators located at the root, 
midspan, and tip. An examination of these records indicated that there 
was a negligible difference in deflections at the root, midspan, and tip 
of the control surfaces for the flight conditions investigated. Therefore 
only the midspan deflection records were used in the analysis. Wing accel 
erations were measured with Statham linear accelerometers. These quanti- 
ties were recorded on Weston 12-channel and Consolidated 18-channel oscil- 
lographs. Self-recording turn meters were used to measure roll and yaw 
rates at the center of gravity and tail. Sideslip angles were measured 
with a vane- type pickup. The instruments were aligned with respect to 
the body reference line and locations of the pertinent instruments are 
shown in figure 2. 


TEST PROCEDURE 


Responses to both aileron and rudder pulses were obtained for a Mach 
number range from 0.6 to 0.82 at an altitude of 35,000 feet. This Mach 
number range is that in which, from an examination of wind-tunnel data, 
the aerodynamic derivatives could be expected to be constant for a given 
flight condition, and thus linear analysis techniques could be used. The 
center-of-gravity location was approximately 21 percent of the mean aero- 
dynamic chord and the gross weight was approximately 115,000 pounds. The 
maneuvers consisted of the pilot applying a sharp aileron or rudder pulse 
with the other controls held fixed, and then allowing the airplane to 
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oscillate for about 25 seconds (yaw damper off). The transient motion 
did not completely subside in this time because of the lightly damped 
Dutch roll mode. Typical time histories are shown in figures 3 and 4. 

In order to show the structural vibrations more clearly, only the first 
portion of several of the records is presented. In addition, complete 
records are shown for roll rate, yaw rate, and side slip- angle responses 
to the rudder pulse to illustrate the motion due to the Dutch roll mode. 
For some of the maneuvers , a small amount of Dutch roll motion was present 
when the control input was applied. This effect has been subtracted out 
of the time histories shown in order to indicate more clearly the response 
to the control-surface motions. 


PREDICTED DYNAMIC RESPONSES 


The dynamic -response characteristics of the airplane are more easily 
interpreted if the time histories are converted to frequency- response or 
transfer-function form. In order to establish the relationships between 
these forms of the dynamic response, a knowledge of the equations of 
motion is necessary. These equations also provide a basis for deter- 
mining approximate relations for transfer-function coefficients and 
aerodynamic derivatives. 


Equations of Motion 


In order to define completely the motions of a flexible airplane, a 
prohibitively large number of degrees of freedom would be needed. Hence, 
for practical purposes, only the most significant degrees of freedom can 
be included. A simplification of particular interest from the standpoint 
of control studies is the use of equations of motion for a rigid airplane 
with the coefficients modified to include the pseudostatic effects of 
flexibility. In the pseudostatic method, as discussed in references 1 
and 2, the structural deformations are assumed to be in phase with the 
aerodynamic and inertial loads resulting from airplane motions. The 
aerodynamic forces arising from these deformations are included by modi- 
fying the derivatives in the airplane equations of motion. 

Three-degree-of-freedom lateral-directional equations of motion with 
respect to stability axes are given in Appendix A. For the flexible air- 
plane appropriate derivatives corresponding to structural deflections due 
to inertial and dead-weight forces resulting from airplane motions have 
been added to the equations. These additional terms are then combined 
with the corresponding rigid airplane inertias to form an effective iner- 
tia for the flexible airplane. As explained in Appendix A, it was found 
to be permissible to neglect a number of the inertial derivatives since the 
cumulative effect of these quantities on the transfer-function coefficients 
was found to be negligible. 
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Aerodynamic Derivatives 


The stability derivatives were obtained by available theory and wind- 
tunnel data. Effects of structural deflections on the derivatives were 
obtained through use of aerodynamic and structural influence coefficients 
in the manner described in reference 1. The chordwise centers of pressure 
of all aerodynamic loads , except loads due to the rudder and ailerons , were 
assumed to be at the quarter-chord line of their respective surfaces. The 
average chordwise center of pressure for loads due to the rudder was esti- 
mated to be the 0.^+8 chord from the pressure-distribution data given in 
reference 5* Although no pressure-distribution data were available for 
the specific wing-aileron plan form, the center of pressure for the aileron 
loads was estimated to be at the O.58 chord from an examination of pressure 
distributions on similar swept-wing plan forms. Fuselage influence coeffi- 
cients used in the analysis were obtained from data given in references 6 
and 7 and are listed in table II together with the wing influence coeffi- 
cients which were obtained from reference 8. Wing stations and equivalent 
weights used in the analysis are shown in figure 5* An example of the 
application of the method to obtain the effect of wing deflections on the 
derivative C 7 is given in Appendix B. Predicted aerodynamic derivatives 

for both the flexible and rigid conditions are listed in table III. 


Transfer-Function Coefficients 


Transfer functions obtained by solving the three-degree-of-freedom 
equations are presented in Appendix A. As noted in this Appendix, some 
of the transfer functions may be simplified for the frequency range of 
interest. After the spiral mode is neglected, rolling response to aileron 
becomes 

_P_ = AafD 2 + aj? + a g ) 

S a (D + D r )(D 2 + CjP + c 2 ) 


The transfer function may also be written as 

p Kp/aJ 1 + ^a p /“a * (P/V 2 ! 

6 a (1 + tD) [1 + 2?D/wn + (D/c^) 2 ] 


For yawing and sideslipping response to rudder motion , both spiral and 
rolling modes may be neglected so that 



10 


NACA TN 387^ 


j, HpD 

8 r D 2 + CjD + c 2 


( 3 ) 


JL = J P 

5 r D 2 + CjD + c 2 


(M 


These quantities may also he expressed as 


r VSr D 

5 r 1 + 2p/u n + (D/c^) 2 


( 5 ) 


J_ = K p/5 r 

5 r 1 + 2^/^ + (D/u^) 2 


The predicted frequency responses may be obtained from the transfer 
functions after iw is substituted for the operator D. The conversion 
of numerator transfer-function coefficients from stability to body axes 
for comparison with frequency responses measured with respect to body 
axes is given in Appendix A. r 


Approximate Equations to Obtain Aerodynamic Derivatives 
From Transfer-Function Coefficients 


The equations for the aerodynamic derivatives were obtained as 
follows: The coefficients C 3 , C 2 , and C-l of the characteristic equa- 

tion given in Appendix A were expressed in terms of their principal aero- 
dynamic derivatives . The coefficients were also expressed in terms of 
the factors of the characteristic equation D r , c x , and c 2 . By equating 
the coefficients of like powers of D, three simultaneous equations were 
obtained with three unknowns, Lp, N r + , and * . A cubic equation 

for Lp was then obtained by eliminating the other unknowns . The deriv- 
ative Lp was found by using an approximate solution for a cubic equation 
together with an approximate equation for the amplitude of the roll to yaw 
ratio of the Dutch roll mode, |£|. Equations for the two other derivatives 
could then be obtained. The equations are 
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Lp = -Dr “ M 
+ = M- " c i 

Np' = c 2 + n(D r - c x ) 


(7) 

( 8 ) 
(9) 


where 


M- = 


1 1 1 (Np - Yep - W.) 


7 1 + D r 2 / C 2 


( 10 ) 


Note that these equations correspond to those given in reference 4 except 
for the quantity p. The additional term p could he considered as a 
measure of the coupling between the Dutch roll and rolling modes. The 
equation used for L^' was obtained from reference 4 


v - Ilk '-/ 1 - v 2 /®* < u > 

Approximate relations for control effectiveness were obtained by substi- 
tuting typical numerical values in the equations given in Appendix A, and 
the following expressions ( analogous to those given in ref. 4) were found 
to be valid for the present airplane. The rudder effectiveness derivative 
was obtained from 


N S r = H r or “Jp 


( 12 ) 


while aileron effectiveness was obtained from 

L S„ = A 3 


(13) 


Aerodynamic derivatives may be found from the quantities on the left-hand 
side of equations (7) through ( 13 ) through use of the definitions given in 
the notation. 


Thus the aerodynamic derivatives may be obtained from transfer- 
function coefficients matched from experimental data, provided that an 

estimate is made of the product of inertia, the derivative C n , and also 

P 

the effective moments of inertia for the flexible airplane. Estimated 
variations in effective inertia characteristics are presented in figure 6. 
The principal effect of flexibility is to reduce the rolling moment of 
inertia. This effect is largely due to the mass of the wing-mounted 
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nacelles. In order to demonstrate the adequacy of the approximate 
equations , predicted transfer-function coefficients were first calculated 
from the predicted derivatives (table III) using the complete expressions 
for the transfer-function coefficients. Then approximate values of the 
derivatives were calculated from the predicted transfer-function coeffi- 
cients using equations (7) to (ll) , and also using the equations with 
^ = 0, which correspond to the equations given in reference 4. Comparisons 
of these derivatives are shown in figure 7 . Also shown are derivatives 
calculated from equations given in reference 4. In all cases the values 
of derivatives using the approximate equations are closer to the estimated 
values than those obtained using the equations from reference 4. This is 
particularly true for the derivative, C nr + . Although not shown, 

a similar comparison was made for the rigid case, and the inclusion of the 
quantity p also improved the agreement with the predicted derivatives. 


REDUCTION OF EXPERIMENTAL DATA 
Frequency Response 


The Fourier integral was used to convert the transient responses to 
frequency- response form. The calculation of the Fourier integral con- 
sisted first of a numerical integration over the first part of the record- 
using a parabolic fit of the time history similar to that described in 
reference 9 . For the control inputs the time histories showed a constant 
value after a certain time and hence the remaining portion of the integral 
was evaluated using an analytical expression for the small step displace- 
ment. The response time histories did not reach a constant value since 
the Dutch roll oscillation was very lightly damped, and analytical expres- 
sions for both a step and a damped sinusoidal oscillation (e.g., ref. 10 ) 
were used to evaluate the remaining part of the integrals. These calcu- 
lations were performed on a card-programmed IBM 650 digital computer. 

All traces were read at 0.05-second intervals except sideslip angle, 
which was read at 0.1-second intervals. In some cases, Fourier transforms 
were obtained from roll and yaw angular- acceleration traces in order to 
define better the corresponding angular rate traces at higher frequencies. 
Results from sideslip records at the higher frequencies are not presented 
because of the lack of well-defined trends. This is probably due to the 
relative attenuation of the higher frequency amplitudes measured by an 
angular displacement instrument as compared with those measured by angular 
rate instruments . Phase-angle corrections were made for the dynamic 
response of the instruments . 
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Transfer-Function Coefficients and 
Aerodynamic Derivatives 


Transfer-function coefficients of the simplified transfer functions 
for p/6 a , r/6 r , and p/& r described in a previous section were obtained 
by matching the response time histories (e.g., figs. 3 and 4 ) by means of 
an analog computer. They were then converted from body axes to stability 
axes. In some cases, small residual oscillations, which were principally 
due to the lightly damped Dutch roll mode, were present when the airplane 
control pulse was applied. Initial conditions were applied to the simu- 
lator to include this motion, although in all cases this had only a small 
effect on transfer-function coefficients. In matching the time histories 
of p/& a on the analog computer, it was found convenient to represent the 
third-order transfer function by parallel networks of first-order and 
second-order transfer functions. This representation is equivalent to 
separating the third-order transfer function into the partial fraction 
form shown in Appendix A. In this way contributions of the Dutch roll 
and rolling modes to the resultant motions could be matched separately. 

Aerodynamic derivatives were then- obtained from the appropriate 
transfer- function coefficients using the approximate equations (eqs. (7) 
to (13)) given in the previous section. 


RESULTS AND DISCUSSION 
Frequency Responses 


In the discussion that follows, effects of the relatively low 
frequency lateral-directional modes will be investigated first. Com- 
parisons will be made between measured and predicted responses to deter- 
mine the frequency range for which the predictions are adequate. Then 
the higher frequency structural modes will be discussed and compared with 
values obtained from the ground oscillation tests of reference 11 . The 
frequency- response data to be presented include: measured frequency 

• responses to aileron and rudder inputs over a fairly wide Mach number 
range, (figs. 8 and 9) > typical comparisons of measured responses at the 
center of gravity and tail (figs. 10 and 11), and typical comparisons of 
measured and estimated frequency responses (figs. 12 and 13). 

Lateral-directional modes .- The ailerons excite principally the 
rolling mode as well as some Dutch roll mode as shown in the roll- rate 
response (fig. 8(a)), which is seen to be of the same form as that indi- 
cated in equation ( 2 ). The rolling mode (a first-order teira) causes the 
net shift of 90° in the phase angle between the frequencies of 0.5 to 5 
radians per second and the slope of the amplitude curve in the frequency 
range of 2 to 7 radians per second, while the peak in the amplitude curve 
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that occurs at frequencies from about 1 to 1.5 radians per second is due 
to the Dutch roll mode (a second-order term) . Responses in sideslip and 
yaw rate are quite small (figs. 8(b) and 8 (c) ) . The rudder mainly excites 
the Dutch roll mode (fig. 9). For both inputs, differences in tail and 
center-of -gravity roll and yaw rates are small in this frequency range. 
This may also be seen from the time histories (figs. 3 and 4). 


The predicted results for the flexible airplane in general agree 
reasonably well with experimental results for frequencies up to about 3 
to 8 radians per second, with the possible exception of wing- tip accel- 
eration (figs. 12(e) and 13(f)). Also some discrepancies are apparent 
in the lowest frequency range (near 0.5 radian/sec) which is the region 
in which the experimental, frequency- response accuracy obtained from pulse- 
type inputs becomes relatively poor, as was discussed in references 1 
and 4. For the aileron responses, the predicted gain for the flexible 
case is much lower than that for the rigid, reflecting the reduction in 
aileron effectiveness. However, the p/6 r responses (fig. 13 ( a )) indi- 
cate a higher gain for the flexible case. This is largely due to the 


higher value of the ratio of 


to the effective moment of inertia in 


the flexible case. Flexibility also reduces the frequency and damping of 
the Dutch roll mode. A discrepancy between theory and experiment in the 
curves for r/6 a (fig. 12 (b)) is that the experimental phase angles 
shift l80° while the predicted values approach a 360° shift. One possi- 
ble explanation is that the smaller phase shift could have been obtained 
from the predicted curve if a more negative value of C n had been 

°a 

assumed. A structural mode could also affect the phase angle in this 
region. 


Structural modes.- A peak occurs in all the frequency responses to 
aileron (fig. 8). at about 13.4 radians per second which will be called 
the first wing antisymmetric bending mode. This compares with l4.4 
radians per second obtained from ground vibration tests (ref. 11). The 
mode is particularly noticeable in the center-of-gravity roll rate and 
wing- tip acceleration responses. In the time history for the wing opti- 
graph, this mode is seen to predominate. Another point of interest is 
that the amplitude of the yaw rate at the tail is greater than that for 
the center of gravity, as may be seen from both the frequency-response 
and time-history curves. The mode frequency remains essentially constant 
for the flight range investigated although there is a slight increase in 
amplitude with Mach number. Measurements of the time history indicated 
that the damping was also essentially constant over the flight range at 
a value of about 0.025 critical. 


The shape of this mode, obtained from measurements of the free 
oscillation amplitudes of the accelerometer and turn-meter traces result- 
ing from aileron inputs, is shown in figure l4. When the curve was deter- 
mined, small corrections for chordwise displacement of the instruments were 
made where necessary to transfer the accelerometer results to the 17- 
percent and 58-percent chord lines corresponding to the front and rear 
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spars, respectively. In the placement of pickups for an autopilot system, 
it is desirable to locate them at points on the structure that are least 
affected by the structural modes in order to avoid introducing undesired 
signals. For a roll gyro, the optimum location for minimizing the angular 
velocity input of the antisymmetric wing mode would be where the wing mode 
shape has zero slope which, as shown in figure lU, is at t] = O .38 near 
the inboard nacelle, while for a linear accelerometer, the optimum loca- 
tion would be at a nodal point, which in figure lU is at q of about 0.55 
just outboard of the nacelle. The wing accelerometer node line obtained 
from the two span positions shows good agreement with that obtained from 
ground vibration tests. 

In the wing- tip acceleration curve (fig. 8(d)), a dip in the amplitude 
and a shift of approximately l80° in the phase angle relative to the 
center-of-gravity roll rate (fig. 8(a)) occurs in a frequency range some- 
what lower than the wing bending frequency. This motion can be illustrated 
as follows: At low frequencies, the relation between wing bending and roll 

angle is as shown in sketch (a). As the frequency increases, the amount of 



wing bending relative to roll angle becomes greater so that the accelera- 
tion at the tip reverses sign with respect to the roll angle (sketch (b)). 



This trend continues until the bending frequency is reached. 

An examination of figure 8 shows that no additional modes were excited 
to any appreciable extent in the measured frequency range above that of the 
wing first antisymmetric bending mode, although two modes were measured 
from ground vibration tests (ref. 11). 

While the rudder pulses also excite the wing antisymmetric bending 
mode, their principal effect is to excite two slightly higher frequency 
modes which are predominantly fuselage modes as may be seen from the 
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responses of yaw and roll rate at the tail and center of gravity (figs. 9 
and 11 ) . The lower frequency, which is predominantly fuselage side bend- 
ing, is 16.5 radians per second while the upper frequency, which is prin- 
cipally fuselage torsion, is 23 radians per second. These frequencies 
compare with ground vibration test values of 18.0 and 24 -. 5 radians per 
second, respectively (ref. 11 ). The prominance of these modes at the 
tail may also be seen from the transient responses due to the rudder 
input (fig. 4 -). 

The frequency responses for roll and yaw rate at the tail due to a 
rudder pulse (fig. 11) also indicate regions of low amplitude at frequen- 
cies below the structural natural frequencies (about 3 'to 8 radians per 
second) . While phase angle for the tail yaw rate shifts l 80 ° out of phase 
with that of the center of gravity, as was the case for the wing at fre- 
quencies below that of the antisymmetric bending mode, the phase angle 
shifts for the tail and center-of -gravity roll rate are approximately the 
same so that they remain in phase. It is of interest to note that by the 
inclusion of only the pseudostatic effects of sideslip, yawing accelera- 
tion, etc., on the tail deflections, these amplitude and phase-angle trends 
are predicted (figs. 13 (b) and 13(d)). While the direction of the 180 ° 
phase-angle shift does not agree for the tail yaw rate, the resultant shift 
is the same. Since the amplitude ratio is quite low in this region, some 
question does exist as to the direction of the phase-angle shifts for both 
the experimental and predicted values. 


Transfer-Function Coefficients 


Comparisons of experimentally determined transfer-function coeffi- 
cients with predicted values (with respect to stability axes) for both 
the flexible and rigid case based on the simplified transfer functions 
(eqs. ( 2 ), (5), and (6 )) are shown in figure 15. In general, agreement 
between experiment and theory for both numerator and denominator is seen 
to be good, as was the case reported in reference 2 for 25,000 feet; 

For the denominator terms , the values of measured damping ratio are 
an average of 0.015 less than the estimated flexible values. While the 
ratios of measured to estimated values are relatively large, the average 
magnitude is approaching zero at the lower values of Mach number, and 
hence the discrepancy is considered to be small. Also some scatter 
occurred in the time constant t. This parameter was difficult to eval- 
uate because it did not have as important an effect on the transient 
response to a pulse as did the other parameters. Flexibility is seen to 
reduce both the frequency and damping of the Dutch roll mode. These 
reductions are caused primarily by the decrease in vertical- tail effec- 
tiveness due to fuselage and tail flexibility having a larger effect than 
the decrease in effective yawing moment of inertia. Flexibility changes 
the time constant t very little. This term depends principally on the 
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ratio of effective rolling moment of inertia to the aerodynamic damping 
in roll. Since the rolling moment of inertia is due largely to the 
nacelle and wing weights , wing flexibility reduces both the aerodynamic 
and inertia loads by about the same amount. 

For the numerator terms also, good agreement is obtained between 
theory and experiment. The scatter that occurred in the measured values 
of £ a resulted from the small effect of £ a on the time histories 
matched with the analog computer. The comparison of the matched curves 
in figures 3 and 4 with the measured time histories indicates that the . 
transfer functions corresponding to the Dutch roll mode only for the rud- 
der pulses, and the rolling and Dutch roll modes for the aileron pulses 
were sufficient to define the time histories adequately. The greatest 
effect of flexibility is seen to be in the reduction in gain for p/b a . 
This reflects the reduction in aileron effectiveness associated with out- 
board ailerons. The predicted values of £ a and w a are slightly differ- 
ent from £ and u) n and approach these quantities at the higher Mach num- 
bers. These differences depend upon the aerodynamic derivatives that 
effect the quantities q and C n , which become relatively small at the 

^a 

higher Mach numbers and corresponding low angles of attack. These . effects 
are reflected in the time histories by. the amount of excitation of the 
Dutch roll mode and also by a reduction in over- all gain. For the limit- 
ing case .of ■ £ a = £ and u a = effects of. the Dutch roll mode are elim- 
inated and the transfer function (eq. ( 2 )) reduces to a first-order form. 


Aerodynamic Derivatives 


Generally good agreement is obtained between predicted and experi- 
mental values of the aerodynamic derivatives (fig. 16) . Differences 
between predicted aerodynamic derivatives for the flexible and : rigid 
cases are generally somewhat greater than the differences between 
transfer-function coefficients since the latter are partially compensated 
by changes in the effective moments of inertia due to flexibility. In 
general, flexibility tends to reduce the values of the aerodynamic deriv- 
atives with the largest reduction occurring for C 7 and C 7 . However, 

L §a P 

flexibility tends to increase the magnitude of C 7 since the. usual 

•3 ' 

decrease in value of a derivative due' to flexibility is more than compen- 
sated by the increase in dihedral due to wing lift. 


CONCLUSIONS ' 


The dynamic lateral-directional responses to aileron and rudder pulses 
of a large flexible airplane for flight conditions of 0.6 to 0.82 Mach num- 
ber and .an altitude of 35 > 000 feet have been investigated and compared with 
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predicted responses "based on the lateral-directional equations of motion. 
Pseudostatic effects of flexibility were included through use of struc- 
tural and aerodynamic influence coefficients. The following conclusions 
have been made: 

1. The manually applied pulse technique was found to be adequate to 
obtain data suitable for defining frequency responses for the frequency 
range from 1 to 25 radians per second. This range includes the relatively 
low frequency Dutch roll and rolling modes, as well as the higher frequency 
wing first antisymmetric bending, fuselage side-bending, and fuselage 
torsion modes. 

2. Reasonably good agreement was obtained in the frequency range 
below about one half of the lowest structural frequency between measured 
and estimated- frequency responses . 

3. The principal effect of flexibility on predicted transfer-function 
coefficients is the reduction in gain of airplane response, to aileron 
motion, reflecting the reduction in aileron effectiveness, and also a 
reduction in damping and frequency of the Dutch roll mode. 

4. The principal lateral-directional derivatives and control- 
effectiveness derivatives, which were evaluated from measured transfer 
functions through use of approximate equations, indicated good agreement, 
in general, with predicted derivatives when the pseudostatic effects of 
flexibility were included. 

5. In the frequency range above half the lowest structural frequency, 
the effects of structural modes were sufficiently large (even at the cen- 
ter of gravity) for the responses measured, so that predictions of the fre- 
quency responses, which included only the pseudostatic effects of flexi- 
bility, were inadequate. 

6. The frequencies of the three structural modes measured in flight 
were within about 10 percent of the frequencies of the corresponding modes 
measured in ground vibration tests. In addition, the wing node line for 
the wing first antisymmetric bending mode measured in flight showed good 
agreement with that obtained from the ground vibration test. 


Ames Aeronautical Laboratory 

National Advisory Committee for Aeronautics 
Moffett Field, Calif., Oct. 8, 1956 
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APPENDIX A • 

PREDICTED AIRPLANE RESPONSE 
Equations of Motion 


The three lateral-directional equations of motion with respect to 
stability axes given in Appendix C of reference k for a rigid airplane 
will be modified for use for a flexible airplane. For the flexible case, 
additional terms must be added to take account of aerodynamic forces 

resulting from structural deflections due to inertial loads (C 7 ... C 7 .., 

6 Cp 6 \|r 

C n > c v.> 811(1 c v ) 811(1 dead- weight loads (Cy m ) . A more complete evalua- 
\|? r 9 

tion of these additional derivatives would also have included other deriv- 
atives dependent on cp, \f, ijr+p, and cp. However, by the substitution of 
typical numerical values in the transfer-function coefficients, the cumu- 
lative effects of these other derivatives on the transfer-function coef- 
ficients were found to be negligible even though some of them were of 
appreciable size (e.g., C n . and C n . were several times larger than 

t P 

C n (t>/2V) for the largest values of M) . The rolling moment, yawing 

moment, and side-force equations, with flight-path angle assumed zero, 
may be expressed as 




V 2 - % ShC l-f - lo Sl,C lp * + 


(- 


t 


-hzP - % SbC l,. D - <lo SbC Z r 2V V* " % ShC l^ = % SlaC K b (A1 ) 


(-= 


I X?P - ^ Sbc 


•o rip 


I? 11 )''* 




I Z D 2 - q o SbC n „D 2 - q Q SbC nr ^ Dj * - c^SbC^P = q^bC^b 


(A2) 



* 1o SC 



♦ + 



(A3) 
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The remaining inertial deflection derivatives were then combined with the 
mass terms to form effective inertias for the flexible airplane. Thus 


t x f " h - io st,c i $ 

V ^ 

q S 

m F = m - ~T C Yg 
<lo s n 

or m C v 

V ' 1 i 

♦ 

Since the derivatives Cy. and Cy_ each represent aerodynamic forces due 

V 3 

to lateral acceleration, they are of equal magnitude, and hence m-p may 
be expressed in terms of either derivative. The three equations (Al) , 
(A2) , and (A3) can then be written in a more convenient form by dividing 
them by Iy , I^, and m^V, respectively, and introducing new symbols 


(f - y>) cp + (-r Xj y - Itd) t - V V 


(a4) 


(- v 2 ' n P D ) 9 + - v>) 


D - NjD ) t - N p p = N 5 6 


(A5) 


“YqjCp + D\|f + (D - Y R )p = Y*6 


3' 


(a6) 


Note that the term (which equals g/v) remains the same as for the 

rigid case since the derivative C v is due to dead-weight loads that 

Xc P 

are distributed in the same manner as the lateral acceleration loads. 


Transfer Functions 


Transfer- function coefficients for cp/5 ^ f/b, and p/5 can be obtained 
by simultaneous solution of the three equations (A*0 , (A5 ) , and (A6) . 



NACA TN 387^ 


21 


Denominator coefficients . - The denominator of the transfer functions, 
A, may be expressed as 

A = D(C 4 D 4 + C 3 D 3 + C 2 D 2 + + C Q ) 

where 



c 2 = Np' + Yp(Lp' + N r ') + LpN r - NpLp 

C ! = L p N P ' ' W + VV*p - N rLp) 

C Q = Ycp(LpN r - Npl*) 

For the moderate angle- of- attack range considered, r v and r^ are 

■ ajp Zjn 

small quantities and hence C 4 * 1. The denominator can then be expressed 
in factored form as 

A = D(D + D g ) (D + D r )(D 2 + c jD + c 2 ) 

where D s and D r represent the spiral and rolling modes, respectively, 
and c x and c 2 are coefficients that define the damping and frequency 
of the Dutch roll mode. 

Numerator coefficients .- With p = Dcp, and r = D\|r, the numerator 
coefficients of the following transfer functions can be expressed in 
terms of the derivatives. 

p _ D(AgD 3 + A ^ 2 + AjD) 

5 A 

where 



' - I P V + % L r - Vr + Wfi' 


A i L s w p " N 5 L P + Y p^ L 8 W r “ + YgCLr^p “ N r L p) 
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r _ D(B S D 3 + B g D 2 + BJ) + B Q ) 
6 " A 


where 


B.. 


B. 


B. 


B, 


- V 

= L 5 N P - V* - Y P V + W 

= Y p ^ ” L 5 N P^ + Y 5^ L 3 N P " ^kp) 
= Wp - W 6 L P ) 

3 _ D(E 3 D 3 + E g D 2 + EjD + Ep) . 
S “ A 


where 

e 3 = Y s( x ' r x F r z F ) 

E 2 = "VV + Lp') - V 

E x = YcpLg* + %% " L s N p + Y 5(Lp N r _ E p^r) 

E o = W* - L g w r) 

Simplifications of the transfer functions .- When transfer functions 
are evaluated from measured data, it is desirable to use as simple a form 
as possible which will still adequately fit the data. Calculations indi- 
cated that the spiral mode factor D s was very small and could be neg- 
lected for the frequency range of interest. Thus the rolling response to 
aileron was simplified to 

p _ A^D 2 + aji) + a-g,) 

5 a (D + D r ) (D 2 + CjD + c 2 ) 


where 
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To determine estimated values for use with the curve fitting of the 
measured responses on the analog computer, a partial fraction form of 
the transfer function was advantageous to use. 


_P_ 

6 a 


Gp 

D + Dp 


+ 


HpD + J 
D 2 + c-jD + c 2 


In obtaining approximate relations for r/b r and p/ 5 r , it was 
desirable first to write the transfer functions in partial fraction 
form. 

y. F Gj* Hj*D + Jj* 

= + + r 

8 r D + D s D + D r D 2 + CpD + c 2 


_p_ 

8 r 


F, 


HpD + Jp 


D + D s D + D r 


D 2 + c jD + C2 


By the substitution of typical numerical values, all numerator terms were 
found to be negligible except H r and Jp. Thus, the transfer-function 
coefficients could be simplified to 


2* H r D 

& r D 2 + c x D + c 2 

J P 

D 2 + c + c 2 
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Transfer functions for the rear part of the fuselage and the wing 
tip .- In addition to the transfer-function coefficients for quantities 
at the center of gravity, the coefficients for r^/S, p^/6 , and n^/8 are 

needed. The evaluation of these quantities includes responses at the 
center of gravity and also responses due to pseudostatic structural 
deflections at the particular location. The equation for yaw rate at 
the rear part of the fuselage is 

ft = £ DA ^t 

• 5 6 8 

where At represents the total change in angle of yaw in radians at the 
t 

rear part of the fuselage due to structural deformations resulting from 
a control input. The following quantities were included in determining 
Ai|r t /5 

A *t _ A ^t P t r ^ A t t jj | A ^t 
5 35 r 8 t 8 8 


For instance A^ /3 represents the pseudostatic change in angle of yaw 
"t 

at the rear part of the fuselage due to structural deformations resulting 
from a unit change in p, 


Similarly, the equation for roll rate at the rear part of the fuselage 


h _ p 
8 8 + 8 


where Acp represents the total change in angle of roll at the rear part 
t 

of the fuselage due to structural deformations resulting from a control 
input and may be expressed as 


^ = + \r + \± + \ 

8 3 8 r 8 t 8 8 
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The equation for acceleration at the left wing tip is 

fVt = _ Dp 

6 8 + 125 

Calculations indicated that the effect of D 2 z /5 was small for the 

wt' 

frequency range up to the first antisymmetric bending frequency and, 
hence, was neglected. 

Conversion of transfer-function coefficients from stability axes to 
body axes .- In order to compare predicted results with measured frequency 
responses which were obtained with respect to body axes, the predicted 
transfer functions were converted from stability-axis to body-axis form. 
Since, for the angle- of -attack range considered, the quantity C 4 was 
essentially unity, only the numerator terms of the transfer function need 
be modified. The equations for conversion from stability axes to body 
axes are (e.g., ref. 


B 


= p cos a - r sin a 


r-n = r cos a + p sin a 


tan (3 b 


tan p 
cos a 


For small angles of attack, the equations may be simplified to 


P B = p - ra 
r B = r + pa 
P B = P 


These relations will also hold for derivatives of the angles. Hence, the 
numerator coefficients for p/5 and r/S may be converted as follows: 


A_- = A.? - B n -a 


b i b = B i + A i a 
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APPENDIX B 


EXAMPLE CALCULATION OF WING CONTRIBUTION TO C, USING 
AERODYNAMIC AND STRUCTURAL INFLUENCE COEFFICIENTS 


A general description of the method used as applied to symmetrical 
loadings was given in reference 1. An example application to obtain 
antisymmetric loadings will be presented here. The equations to be used 
will be given first, and then a numerical example to obtain C 1 will 
be shown. P 

The rolling-moment coefficient resulting from an arbitrary antisym- 
metric angle-of-attack distribution may be expressed as a power series 
of the dynamic pressure q. 


— Ct + C-7 q + C 
l R l A 


* C l/ + 


(Bl) 


where 


C ^ rolling-moment coefficient for the flexible wing 


C 7 rolling-moment coefficient for the rigid wing 

l R 


Cj q increment in rolling-moment coefficient resulting from structural 
A deflections due to the Cj^ loading 


Cj I 2 
t B 


increment in rolling-moment coefficient resulting from structural 
deflections due to the C j loading 


For the wing structure considered in the present example, the ratios 
of all terms after the first one, -C^ /C^ , -C^ /C^ , etc., are essen- 

A B B C 

tially equal to a single constant k. This result may be interpreted to 
mean that the successive wing-deflection shapes that produce these loads 
are essentially the same, and, hence, this portion of the wing loading 
may be represented as a single-degree-of-freedom system. Thus the incre- 
mental rolling-moment coefficient due to this portion of the wing deflec- 
tions may be represented by the expression C^q/l+kq and it will be 

valid for large positive values of the quantity kq even though the series 
would not converge for positive kq greater than 1. The equation for 
rolling-moment coefficient for the flexible wing will then be 
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C 


l 


F 



1+ ^VV 

1 + kq 


(B2) 


The quantities in equation (B2) will now be determined through use of 
aerodynamic and structural influence coefficients. 

Antisymmetric aerodynamic influence coefficients were obtained through 
use of reference 12. First, angle-of -attack influence coefficients, 

which represent the angle of attack at station v due to a continuous 
loading function having a unit value at station n and zero value at the 
other stations, were obtained from reference 12. These were then con- 
verted to loading influence coefficients A yn (i.e., the loading c^c/2b, 

at station, n, due to a continuous angle -of- attack distribution with a 
unit value at station v and zero values at the other stations) by solving 
for the span loadings for a unit angle of attack at one control station. 
This was done for each of the spanwise control stations 1, 2, and 3 located 
at T) = 0.924, 0.707 y and O.383, respectively. The resulting rolling- 
moment coefficient for a unit angle of attack at a control station, v, can 
then be calculated using the following equation, which has the form of 
equation 15 in reference 12. 


o Iv = [a v2 + 0.707 (a v1 + a v3 ) ] 


(B3) 


The wing structural influence coefficients are listed in table II. 
For use with the aerodynamic influence coefficients, they were converted 
to the form of a change in streamwise angle of attack in radians, due to 
a 1000- pound load at a given percent chord position. All aerodynamic 
loads were assumed to act along the quarter- chord line except for loads 
due to the ailerons. The change in streamwise angle of attack at sta- 
tion m due to the loading resulting from an angle of attack at station 
v, with q =' 1 psi, may be expressed as 

■ 3 

(€o) vm = ^ 1000 X A vnWn • W 

n=i 


with 


m = 1, 2, 3 

v = 1, 2, 3 
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where 

A vn aerodynamic influence coefficients for the loading function at 
station n due to an angle of attack at station v 

Sjmj structural influence coefficients expressing the change in angle 
of attack at station m in radians due to a 1000-pound load 
at the quarter- chord station n 

I n integrating factors with the following values 

I x = 0.1502 

1 2 = 0.2776 

1 3 = O.3628 

The aerodynamic loading due to (e n ) is then obtained from the summation 

u vm 

3 

V n = X ^ € o^vnAnn ( B 5) 

m=i 

from which the incremental rolling-moment coefficient may be obtained. 

( C 0 V ■ 16 { (Aa> V2 + 0 -T 07 [(A A ) V1 + (A A ) J } (B6) 


The next incremental rolling-moment coefficient 

resulting from the loading coefficients (A^) 
manner . 



due to deflection 


is calculated in a similar 


Thus, the rolling-moment coefficients resulting from the initial angle 
of attack, and from the first and second twist distributions can be calcu- 
lated from the following equations 
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c u ■ L (0,“- 

V=1 

3 

% = X ( c 0 v av 

V=1 v 


(B8) 


(B9) 


The rolling-moment coefficient for the flexible wing for any desired 
angle-of -attack distribution, a v , is then obtained by substituting values 
from the above equations into equation (B2). The reference rolling-moment 
coefficients with the aerodynamic loads at the quarter chord were calcu- 
lated through use of the previous equations (with q = 1 psi) and are 
tabulated in the following table: 


V 

1 

2 

3 

c l 
L v 

0.120 

0.298 

0.326 

C 7 

1 Av 

-.0376 

-.0698 

-.0376 

% 

B v 

.OO856 

.01553 

.00796 


In order to determine the rolling-moment coefficient for a specific 
derivative, the wing angle -of -attack distribution must be known. For a 
unit sideslip angle, the wing contribution to may be separated into 

loadings due to two angle -of- attack distributions (ref. 13): (1) a con- 

stant spanwise angle-of-attack distribution due to the difference in 
effective velocity acting on each wing panel for the yawed attitude and 
(2) an angle-of-attack distribution proportional to the wing dihedral 
angle, which for this case may be considered to be due entirely to the 
upward deflection of the wings resulting from the level flight symmetrical 
air loads. 


For the rigid airplane, the contribution to due to the differ- 


ential velocities acting on each wing panel was estimated from reference lit- 
as 

equations (B7), (B8) , (B9 ) , and (B2) gives 


C Zp /C L = -0.129. Substituting values for a unit angle of attack into 
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'Pi 


= -0.129 1 


0.1451 
0.744 q 
, 0.0321 ; 
+ 0.1451 ^ 


= -0.129 1 - 


0.195q 

1 + 0.221q, 


The dihedral angle , measured in a plane parallel with the YZ plane , 
was calculated through use of the influence coefficients for symmetrical 
loads described in reference 1 and, for the flight range of interest, could 
be expressed as 


Vi 


C L 


0.l035q. 


The calculated distribution of dihedral angle, normalized with respect to 
the value at the tip, is given in the following table: 


4 

0.383 

0.707 

0.924 

r 

O 

OJ 

ir\ 

• 

.900 

.992 


For the wing at a unit angle of sideslip, these values of dihedral angle 
represent changes in the wing angle of attack, and the rolling-moment 
coefficient due to dihedral angle may be determined from equations given 
previously. 



-(0.10350) (0.557) 



0.1197 

0.557 

, 0.02664 

1 + 0 

0.1197 y 


= -0.05770 



0.2l4q \ 

1 + 0.222q/ 


Since only first-order effects have been considered, the changes in 
dihedral angle due to antisymmetric loads have been neglected. 
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The total value for the wing contribution to C, becomes 

. 4 P 




-0.129 




__0 1 195i_A 

1 + 0.221q ) 


0.05771 



0.2lAq \ 
1 + 0.222q/ 


(BIO) 


While the two load distributions yield somewhat different values of the 
ratio c Z^/ ( -'Zp^ the va l ues °f k are essentially the same. This was 


also found to be true for other types of load distributions. Note that 


the values of C-, 


obtained from equation (BIO) are not the same as 


those given in table III since only the contribution of the wing has 
been considered here. 
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TABLE I.- PHYSICAL CHARACTERISTICS OF THE TEST AIRPLANE 


Wing 

Area, sq ft 1428 

Span, ft 

Aspect ratio ■ • 9*^3 

Taper ratio 0.42 

Mean aerodynamic chord, ft ' . 13-0 

Sweepback of quarter- chord line, deg 35 


Airfoil thickness ratio (parallel to body center line). 


percent 

Dihedral, deg 0 

Ailerons 

Area aft of hinge center line (each), sq ft 53.8 

Aileron span to wing span ratio 0.405 

Average aileron chord to wing chord ratio 0.26 

Vertical tail 

Area (including dorsal), sq ft 227 

Span, ft 18.7 

Aspect ratio 1.5^ 

Taper ratio 0.34 

Mean aerodynamic chord, ft 13.02 

Sweepback of quarter- chord line , deg 3^.9 

Airfoil thickness ratio (parallel to body center line), 

percent . 10 

Distance l/4-chord wing M.A.C. to l/4-chord vertical tail 

M.A.C., ft ^6-5 

Rudder 

Area aft of hinge center line, sq ft 51*2 

Average rudder chord to vertical tail chord ratio 0.30 

Average gross weight, lb 115,000 

Moment of inertia about the X axis, slug-ft 2 1,074,000 

Moment of inertia about the Z axis, slug-ft 2 2,306,000 


Inclination of principal longitudinal axis with respect 
to fuselage reference line, deg 


- 2.6 



TABLE II.- STRUCTURAL INFLUENCE COEFFICIENTS 
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ite: Wing station locations are shown in figure 5» The letter F denotes 

front spar location; R denotes rear spar location. 
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TABLE III.- LATERAL- DIRECTIONAL STABILITY DERIVATIVES PREDICTED FROM WIND- 
TUNNEL DATA AND AVAILABLE THEORY; h = 35,000 FEET, W = 115,000 POUNDS 


Quantity 

Flexible 

Rigid 

M 

0.5 

0.6 

0.7 

0.8 

0.5 

0.6 

0.7 

0.8 


87.2 

125.5 

171.0 

223.1 

87.2 

125.5 

171.0 


C L 

.922 

.642 

.471 

.361 

.922 

.642 

.471 

.361 

C 7 

Ip 

-.430 

->29 

-.427 

-.426 

-.490 

-.507 

-.540 

-.570 

On 

n P 

- . 0960 

- .0660 

-.0470 

-.0350 

-.099^ 

- . 0692 

-.0507 

-.0389 

C 7 

L r 

• 215 

.136 

.090 

.060 

.264 

.184 

• 135 

.104 

c n 

n r 

-.142 

-.137 

-.13^ 

-.130 

-.150 

-.150 

-.150 

-.150 

% 

-.170 

- ,l4l 

-.122 

-.112 

-.157 

-.129 

-.113 

-.101 

% 

.110 

.106 

.102 

.098 

.120 

.120 

n 

.120 

Cy 

Y P 

-.575 

-.566 

-.555 

-.544 

-.600 

-.600 

-.600 

-.600 

c 7 

.0052 

.0090 

.0109 

.0122 

.0057 

.0095 

.0119 

.0135 

C n fi 

S r 

-.0724 

-.0700 

-.0672 

-.0645 

-.0778 

-.0778 

-.0778 

-.0778 

C Y 

Y S r 

.183 

.180 

.175 

. .171 

.194 

.194 

.194 

.194 

C 7 

.059^ 

.0528 

.0462 

.0392 

.0770 

.0770 

.0770 

.0770 

c n 

-.0103 

- . 0064 

-.0041 

-.0027 

-.0133 

-.0093 

-.0068 

-.0052 

'x 

-.116 

-.059 

-.005 

.031 

-.103 

-.047 

-.010 

.014 

T Jl 1 

-.048 

-.021 

-.005 

.007 

-.047 

-.020 
1 

-.005 

.007 
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Figure 1.- Photograph of the test airplane. 
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Figure 2.- Two- view drawing of test airplane. 
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t, sec 

(b) Time histories of 0 , r cg , r^ , n v t> an( ^ z wt* 

B ® 

Figure 3.“ Concluded. 



, deg p t , radians/ sec p rn , radians/sec 
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— Results from analog computer 




t, sec 


(a) Time histories of 6 r , p , p, , and 3 . 

C °B X B 

Figure 4.- Typical time histories for rudder pulse; M = 0.71. 





, radians / sec r , radians / sec 
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Figure 5 .- Wing stations and weight distribution used in analysis. 
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50 .55 .60 .65 .70 .75 .80 .85 

M 

Figure 6.- Influence of flexibility on effective inertia characteristics. 
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Predicted derivatives 

Derivatives calculated by equations from this report, 
Derivatives calculated by equations from reference 4 


Figure 7.- Comparison of predicted aerodynamic derivatives for the 
flexible airplane with those calculated from predicted transfer 
function coefficients using approximate equations . 















Phase angle, <J> rcg , deg Amplitude ratio, -y-2. t per sec 
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.5 12 5 10 20 30 

Frequency, cj, radians/sec 

(b) Yawing velocity at the center of gravity. 

Figure 8.- Continued. 





Phase angle, <J>£ , deg Amplitude ratio, 





Phase angle, <£> n , deg Amplitude ratio, g/deg 


50 NACA TN 387^ ■ 




Frequency, w, radians/sec 
(d) Wing-tip acceleration. 
Figure 8 .- Concluded. 
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Figure 9 *~ Frequenc; 


rniSKi 
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.5 12 5 10 20 30 


Frequency, a>, radians/sec 
(c) Yawing velocity at the tail. 
Figure 9 .- Continued. 
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Frequency, cj, radians/sec 
(d) Sideslip angle. 


Figure 9 -“ Continued. 




Phase angle, <f> n , deg Amplitude ratio, -j- , g/deg 
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Frequency, w, radians/sec 
(e) Wing-tip acceleration. 
Figure 9 .- Concluded. 
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Frequency, cu, radians/sec 


Figure 10. _ Comparison of the measured, yawing velocity frequency responses 
at the center of gravity and the tail to an aileron input; M = 0.71. 




Phase angle, <$ p , deg Amplitude ratio, -y- , per sec 
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(a) Rolling velocity. 

Figure 11.- Comparison of measured frequency responses at the center of 
gravity and the tail to a rudder input; M = 0 . 71 . 
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Frequency, w, radians/sec 


(a) Rolling velocity at the center of gravity. 

Figure 12.- Comparison of experimental and predicted frequency responses 

to aileron input; M = 0.71. 




60 


MCA TN 387^ 




.5 12 5 10 20 30 

Frequency, co, radians/sec 

(b) Yawing velocity at the center of gravity. 

Figure 12.- Continued. 









Phase angle, 
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(d) Sideslip angle. 
Figure 12.- Continued. 




Phase angle, ® n , deg Amplitude ratio, — , g/deg 
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Frequency, ai, radians/sec 
(e) Wing- tip acceleration. 
Figure 12.- Concluded. 





Phase angle, <£ p , deg Amplitude ratio, |-j— |, per sec 
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(a) Rolling velocity at the center of gravity. 

Figure 13.- Comparison of experimental and predicted frequency responses 

to rudder input; M = 0.71. 
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13.- Continued 
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Experiment 
Predicted flexible 
Predicted rigid 
Results from analog 
computer matching 
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Frequency, cu, radians/sec 

(c) Yawing velocity at the center of gravity. 
Figure 13 .- Continued. 
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(d) Yawing velocity at the tail. 


Figure 13 .- Continued 




Phase angle, <£„ deg Amplitude ratio, 
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Frequency, a>, radians/sec 

(e) Sideslip angle. 
Figure 13 .- Continued. 







Node line 



Figure l4 . - First antisymmetric wing bending free oscillation mode shape; 
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